# The website says the Amcerican Community Survey (ACS) has their most recent data set as 2014-2018. The ACS includesvariables on social, economic, housing, and demographic aspects. It has 25,070 variables.
var <- load_variables(year = 2018, dataset = "acs5")# Finding the demographic variables and renaming them in a list
demo_varsacs5 <- c(
black_population = "C02003_004",
white_population = "C02003_003",
nativeam = "B02010_001",
asian = "B02011_001",
hislat = "B03002_001",
hawpac = "B02012_001",
poverty = "B17001_002",
veteran = "B21001_002"
)Census tracts are small, relatively permanent statistical subdivisions of counties designed to be relatively homogeneous units with respect to population characteristics, economic status, and living conditions. They are updated by local participants prior to each decennial census as part of the Census Bureau’s Participant Statistical Areas Program. Census tracts generally have a population size between 1,200 and 8,000 people, with an optimum size of 4,000 people. Census tract boundaries never cross county lines
# We have to use tracts so that we can get geographical information using sf are leaflet
nj_demo_tracts <- get_acs(
state = "NJ",
geography = "tract",
variables = demo_varsacs5,
# Use 2018 as year for more recent data
year = 2018,
geometry = TRUE,
keep_geo_vars = TRUE
) %>%
# 4326 is most common map form
st_transform(4326)# Function to create tracts
tracts_func <- function(data, county_n, var_n) {
data %>%
# Filters the county
filter(COUNTYFP == county_n) %>%
# Filters for the variable from our variable list
filter(variable == var_n)
}#Function to pull variable from all tracts in TN SD-26
alco_func <- function(var_n){
nj_demo_tracts %>%
filter(COUNTYFP == "003" #bergen
| COUNTYFP == "031" #passaic
| COUNTYFP == "013" #essex
| COUNTYFP == "027" #morris
| COUNTYFP == "037" #sussex
| COUNTYFP == "017" #hudson
| COUNTYFP == "039" #union
| COUNTYFP == "041") %>% #warren
filter(variable == var_n)
}#function to color the map. takes in string color pallette and mapping variable
pal_func <- function(color, var_name){
colorNumeric(
palette = color,
domain = var_name$estimate,
reverse = TRUE
)
}# Function for creating leaflet map
map_func <- function(data, title, title2) {
data %>%
leaflet(width = "100%") %>%
addProviderTiles(provider = "CartoDB.Positron") %>%
addPolygons(
label = ~ str_extract(NAME.y, "^([^,]*)"),
stroke = FALSE,
smoothFactor = 0,
fillOpacity = 0.7,
color = ~ pal(estimate)
) %>%
addLegend("bottomright",
pal = pal,
values = ~estimate,
title = title,
opacity = 1
) %>%
addControl(title2, position = "topright")
}albl_tracts <- alco_func("black_population")# Creating pallette to use with leaflet
pal <- pal_func("viridis", albl_tracts)
map_func(albl_tracts, "Black Population", "Black population Northern NJ")alhs_tracts <- alco_func("hislat")pal <- pal_func("viridis", alhs_tracts)
map_func(alhs_tracts, "Hispanic/Latino Population", "Hispanic/Latino population in Northern NJ")alwh_tracts <- alco_func("white_population")pal <- pal_func("viridis", alwh_tracts)
map_func(alwh_tracts, "White Population", "White population TN SD-26")alas_tracts <- alco_func("asian")pal <- pal_func("viridis", alas_tracts)
map_func(alas_tracts, "Asian Population", "Asian population in Northern NJ")alna_tracts <- alco_func("nativeam")pal <- pal_func("viridis", alna_tracts)
map_func(alna_tracts, "Native American & Alaskan Population", "Native American & Alaskan population in Northern NJ")alha_tracts <- alco_func("hawpac")pal <- pal_func("viridis", alha_tracts)
map_func(alha_tracts, "Hawaiin and Pacific Islander Population", "Hawaiian and Pacific population in Northern NJ")allpov_tracts <- alco_func("poverty")pal <- pal_func("plasma", allpov_tracts)
map_func(allpov_tracts, "Population Living in Poverty", "Number of Impoverished People in Northern NJ")allvet_tracts <- alco_func("veteran")pal <- pal_func("magma", allvet_tracts)
map_func(allvet_tracts, "Veteran Population", "Number of Veterans in Northern NJ")fybl_tracts <- tracts_func(nj_demo_tracts, "003", "black_population")pal <- pal_func("viridis", fybl_tracts)
map_func(fybl_tracts, "Black Population", "Black population of Bergen County")fylat_tracts <- tracts_func(nj_demo_tracts, "003", "hislat")pal <- pal_func("viridis", fylat_tracts)
map_func(fylat_tracts, "Hispanic/Latino Population", "Hispanic/Latino population of Bergen County")fywh_tracts <- tracts_func(nj_demo_tracts, "003", "white_population")pal <- pal_func("viridis", fywh_tracts)
map_func(fywh_tracts, "White Population", "White population of Bergen County")fywh_tracts <- tracts_func(nj_demo_tracts, "003", "asian")pal <- pal_func("viridis", fywh_tracts)
map_func(fywh_tracts, "Asian Population", "Asian population of Bergen County")fypov_tracts <- tracts_func(nj_demo_tracts, "003", "poverty")pal <- pal_func("plasma", fypov_tracts)
map_func(fypov_tracts, "Population Living in Poverty ", "Poverty in Bergen County")fyvet_tracts <- tracts_func(nj_demo_tracts, "003", "veteran")pal <- pal_func("magma", fyvet_tracts)
map_func(fyvet_tracts, "Veteran Population", "Veteran population of Bergen County")hnbl_tracts <- tracts_func(nj_demo_tracts, "031", "black_population")pal <- pal_func("viridis", hnbl_tracts)
map_func(hnbl_tracts, "Black Population", "Black population of Passaic County")hnlat_tracts <- tracts_func(nj_demo_tracts, "031", "hislat")pal <- pal_func("viridis", hnlat_tracts)
map_func(hnlat_tracts, "Hispanic/Latino Population", "Hispanic/Latino population of Passaic County")hnwh_tracts <- tracts_func(nj_demo_tracts, "031", "white_population")pal <- pal_func("viridis", hnwh_tracts)
map_func(hnwh_tracts, "White Population", "White population of Passaic County")hnwh_tracts <- tracts_func(nj_demo_tracts, "031", "asian")pal <- pal_func("viridis", hnwh_tracts)
map_func(hnwh_tracts, "Asian Population", "Asian population of Passaic County")hnpov_tracts <- tracts_func(nj_demo_tracts, "031", "poverty")pal <- pal_func("plasma", hnpov_tracts)
map_func(hnpov_tracts, "Population Living in Poverty", "Poverty in Passaic County")hnvet_tracts <- tracts_func(nj_demo_tracts, "031", "veteran")pal <- pal_func("magma", hnvet_tracts)
map_func(hnvet_tracts, "Veteran Population", "Veteran population of Passaic County")